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Abstract 

Using a sample of 3 x 10® i^i^S) decays collected by the CLEO III and CLEO-c detector con- 
figurations, we present results of a study of Xca ^'^d Xc2 decays into ryry, 7777', and rjvl final states. 
We find S(xcO ^ w) = (0.31 ± 0.05 ±0.04 ±0.02)%, B(xcO ^ wO < 0.05% at the 90% confidence 
level, and i?(xcO ^ W) = (0-17 ± 0.04 ± 0.02 ± 0.01)%. We also present upper limits for the 
decays of Xc2 into these final states. These results give information on the decay mechanism of Xc 
states into pseudoscalars. 

PACS numbers: 13.25.Gv, 14.40.Gx 
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In the standard quark model, the Xcj (J = 0,1,2), mesons are cc states in an L = 1 
configuration. As they cannot be produced directly in e~^e~ collisions, they are less well 
studied than the states. However, 4'{2S) IXcj decays yield many Xcj mesons, and 
e^e~ '0(25') is a very clean environment for Xcj investigation. In this paper, we con- 
centrate on the two-body decays of the Xco and Xc2 into rjr], rji]', and rj'r]' final states^. 
Knowledge from these decay rates leads to information on the quark and gluon nature of 
both the Xc parents and their pseudo-scalar daughters, and a greater understanding of the 
decay mechanisms of the Xc mesons p,]. 

The data presented here were taken by the CLEO detector operating at the Cornell 
Electron Storage Ring with a center of mass energy corresponding to the ipi^S) mass of 
3.686 GeV/c^. The data were taken with two different detector configurations. An integrated 
luminosity of 2.74 pb~^ was collected using the CLEO III detector configuration ^], and 
2.89 pb~^ using the CLEO-c detector configuration [3|. The total number of ?^(25') events 
is calculated as 3.08 x 10^, determined according to the method described in The vital 
detector component for this analysis, the Csl crystal calorimeter is common to the two 
configurations and has an energy resolution of 2.2% at 1 GeV, and 4.0% at 100 MeV. 

We detect our rj candidates in the decay modes 77, tt+tt^tt^, and 7r+7r~7, and our rj' 
candidates in the modes 771+71" and rjir'^TT' , where the t] is reconstructed using the rj decay 
modes previously listed. We look for all combinations of distinct rj and rj' candidates, with 
the exception of the case where both 77 candidates in the event decayed to 77. This last 
case was excluded because its detection depends on the operation of an all-neutral trigger 
whose efficiency varied over the running period because of hardware changes. All our final 
states thus contain at least two charged particles, and the trigger efficiency for these events 
is close to 100% after all other requirements. 

We define photon candidates as clusters in the Csl having a shower shape consistent with 
being due to a photon. We form 17 ^ 77 candidates, from pairs of photon candidates, using 
the event vertex found from the charged tracks as the position of origin of the photons. We 
kinematically constrain the two-photon combination to the known t] mass, and those with a 

< 10 (for one degree of freedom) for this fit are retained for further analysis. The same 
procedure is used for candidates. We then proceed to make r] — > 7r+7r~7r'' candidates by 
constraining the 7r+7r^7r'' combinations to the nominal 1] mass and again requiring a of 
less than 10 for one degree of freedom. Similarly r]' candidates are built, either from the 
?77r"'"7r~ combinations, or from combining 7r"''7r~ and 7 candidates (with >50 MeV), mass 
constraining the resultant 4-momentum to the f]' mass and requiring them to have a < 10 
for the one degree of freedom of this constraint. The mass resolution of the reconstructed 
— > 77 combinations is around 12 MeV/c^ before the mass constraint, and the analogous 
numbers for rj — > 7r+7r^7r'^, 77' 77r+7r~, and rj' r]TT~^TT~ are 3 MeV/c^, 7 MeV/c^, and 2.5 
MeV/c^ respectively. 

If we have two distinct r]^'^ candidates, we combine them into a Xc candidate. At this 
stage of the analysis, the signal to noise ratio is poor, and the invariant mass resolution 
of the Xc is around 15 MeV/c^. We then search for any unused photon in the event, and 
add that to the Xc candidate to make a ip{2S) candidate. This ip{2S) is then kinematically 
constrained to the 4-momentum of the beam, the energy of which is taken as the known 
ip{2S) mass, and the momentum is non-zero only by a tiny amount due to the finite crossing 



We do not consider Xci decays into these final states, as they are forbidden by simple spin-parity conser- 
vation. 
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angle 2 mrad per beam) in CESR. To make our final selection, we require the ipi^S) 
candidate to have a of less than 25 for the 4 degrees of freedom for this fit; this cut rejects 
many background combinations. This kinematic fit greatly improves the mass resolution of 
the Xc candidate. 

To study the efficiency and resolutions, we generated Monte Carlo samples for each Xc 
state into each final state, using a GEANT-based detector simulation [g']. The simulated 
events have a distribution of (1 + A cos^ 6), where 6 is the radiated photon angle relative to 
the positron beam direction, and A = 1 for the Xco and A = 1/13 for the Xc2, in accordance 
with expectations for an El transition. The mass resolution and efficiencies are shown in 
Table I. The resolutions varied a little with the detector configuration and decay chain, 
but are approximated by single Gaussian functions. The efficiencies shown include all the 
relevant branching fractions of the 1] and rj'. 

TABLE I: Efficiencies (in %) and resolutions (in MeV/c^) obtained from analysis of Monte Carlo 
generated events. 



Mode 


XcO 


Xc2 




Efficiency 


Resolution 


Efficiency 


Resolution 




% 


MeV/c2 


% 


MeV/c2 


rjrj 


5.46 


7.6 


5.60 


6.3 


rjrj' 


5.47 


7.0 


5.47 


5.8 


rj'rj' 


4.69 


5.4 


4.67 


4.9 



The final mass plots are shown in Fig. 1. Clear peaks are found for the decays Xco 
and XcO v'v\ but the other four decays under consideration show no significant signals. 
These plots are each fit with two signal shapes comprising Breit-Wigner functions convolved 
with Gaussian resolutions, together with a fiat background term. The masses and widths 
of the Breit-Wigner functions were fixed at their standard values [?!, and the widths of 
the Gaussian resolution functions fixed at the values shown in Table I. We find signals of 
47.8 ± 7.7 events for Xco VV^ and 22.7 ± 5.3 events for Xco v'v'- To find the limits on 
the number of events for each of the other decay modes, we plot the probability function for 
a series of different signal yields, and place the limit so that 90% of the physically allowed 
region is below its value. The signal yields and limits are summarized in Table II. Our 
detected signal events in Xco decays have a 6 angle distribution consistent with the value of 
A = 1 as defined above. They also make a uniform band of events in the M'^{'jrji) versus 
M^{lV2) Dalitz plot, as expected for the decay of a spin-0 particle. 

To convert the yields to branching fractions, we divide by the number of ip{2S) events 
in the data sample multiplied by the detector efficiency and by the branching fractions 
for 4'{2S) into XcJ, for which we use the CLEO measurements of B{%Ij{2S) IXco) = 
9.22 ± 0.11 ± 0.46%, and B{^{2S) 7x^2) = 9.33 ± 0.14 ± 0.61% 

The systematic uncertainties in the branching fractions are summarized in Table III. 
The systematic uncertainties due to fitting the peaks were evaluated by floating, in turn, 
the mass, instrinsic width, and Gaussian resolution width. The maximum change in yield 
for our highest-statistics signal (xco ~^ VV) from this process gives us our measure of this 
systematic uncertainty, 5%. The requirement on the of the constraint to the beam 4- 
momentum has been checked by changing the cut and noting the change in the yield. This 
has been done for the modes under investigation here, and also in other modes such as 
rjn^n^ that have higher statistics. Based on this study we place a systematic uncertainty of 
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3.5% on the efficiency of this requirement. The uncertainties due to track reconstruction are 
small and well understood. The biggest systematic uncertainty is that due to the photon 
reconstruction, which is set at 2% per photon. This was derived from a series of studies of 
photon reconstruction in well understood decays such as ipi^S) J/ipTi^n'^. For the two 
signals Xco ~^ VV Xc2 VV' add the systematic uncertainties in quadrature, except 
for the one due to the ipi^S) —> 'jXc branching fractions which we keep separate, when 
calculating the final branching fractions. For evaluating the limits in the cases where there 
is no significant signal, we take the probability density function, and divide each entry by the 
efficiency smeared by the total systematic uncertainty, and find the branching fraction that 
includes 90% of the total area. A rough check on our ability to find the efficiency of the entire 
decay chain is possible by measuring the number of events consistent with the decay chain 
ip{2S) — >• rjJ/ip, Jip —5- 777'. This has the same final state particles as ^'(25') ^Xc-i Xc ~^ VV'- 
Using the Review of Particle Properties 0] values of the branching fractions, we expect to 
see 19.2 ± 2.4 events in this mode, and we observe 20.8 ± 4.7, which is consistent with this 
expectation. 

TABLE II: Signal yields and branching fraction results for each decay mode. The uncertainties 
are statistical, systematic due to this measurement, and systematic due to the '4>{2S) — > Xcj, 
respectively. The limits on the branching fractions include all systematic uncertainties. 



Mode 


XcO 


Xc2 




Yield 


B.F. (%) 


Yield 


B.F. (%) 


rjf] 


47.8 ± 7.7 


0.31 ±0.05 ±0.04 ±0.02 


< 7.2 


< 0.047 


rjf]' 


< 7.5 


< 0.050 


< 3.4 


< 0.023 


rj'r]' 


22.7 ±5.3 


0.17 ±0.04 ±0.02 ±0.01 


< 4.0 


< 0.031 



TABLE III: Systematic uncertainties expressed in percent. 



Mode 




Trig. 


Fit 


Tracks 


7 


MC stats 


Cuts 


Total 




3 


1 


5 


1.4 


10 


3 


3.5 


12 


7777' 


3 


1 


5 


2.1 


8 


3 


3.5 


11 


ri'i]' 


3 


1 


5 


2.8 


6 


3 


3.5 


10 



Our results can be interpreted in the model of Qiang Zhao [ij], who predicts these decay 
rates as a function of the QCD parameter, r, which is the ratio of doubly- to singly-OZI- 
suppressed decay diagrams. The measurements are consistent with small values of r, which 
indicates that it is the singly suppressed OZI diagram that dominates these decays. We 
also present a limit on the branching fraction for Xc2 decays to rjr] that is tighter than that 
of BES [sl, and the ffist limits for Xc2 decays into r]")]' and t]'-)]'. The limit for Xc2 ~^ VV is 
slightly below the theoretically predicted branching fraction for r = 0, but all three limits 
are consistent with r being small and negative. More definitive conclusions require a larger 
sample of Xc2 decays. 

In summary, we find a branching fraction for Xco VV of (0.31 ± 0.05 ± 0.04 ± 0.02)%. 
This is larger than, but consistent with, two lower-statistics measurements of BES [sj and 
E-835 [9]. We make the ffist measurement of Xco ^ v'v' oi (0.17 ± 0.04 ± 0.02 ± 0.01)%. We 
find no signal for Xco VV' ^'^d are able to set an upper limit on this branching fraction 
of < 0.05% at the 90% confidence level. We also set limits of the branching fraction for Xc2 
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decaying to 7777, rjr]', and ri'r]' of 0.047%, 0.023% and 0.031% respectively. Our results imply 
that singly OZI-suppressed diagrams dominate in these decays. 



FIG. 1: Invariant mass distributions for (a) rjij, (b) ijij' , and (c) rj'r]' . The fits are described in the 
text. 
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